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We report the first observation of coherent back-scattering (CBS) of light in a transverse photonic 
disorder. The CBS peak is recorded in the far-field, at a fixed propagation time set by our crystal 
length, and displays a contrast approaching the ideal value of 1, which proves good coherence of 
transport in our system. We study its dynamics for increasing disorder strength, and find a non- 
monotonous evolution. For weak disorder, the CBS signal increases, and the asymmetry of the 
momentum distribution becomes inverted compared to the initial condition. For stronger disorder, 
we observe a resymmetrization of the momentum distribution, confirmed by numerical simulations, 
and compatible with the onset of strong (Anderson) localization. 

PACS numbers: 


Introduction. Coherent Back Scattering (CBS) is a 
fundamental effect in mesoscopic physics, by which a 
wave of initial momentum ko, multiply scattered elasti¬ 
cally in a disordered medium, has a probability enhanced 
by a factor 2 to end in state —kg, compared to any other 
k within the elastic scattering circle |k| = |ko|. This 
effect is due to constructive interference of the reversed 
pairs of scattering paths, as shown in Fig. [^b, and is inti¬ 
mately connected with weak localization (WL) [T]. It was 
first predicted for electrons [2], and several studies could 
relate observed phenomena in electron transport, such as 
negative magnetoresistance, to WL [3]. However, direct 
observations of the CBS peak were obtained only using 
macroscopic classical waves such as light mm , acous¬ 
tic waves 017 ] , seismic waves [8], and recently ultracold 
atoms [9]. Besides proving the coherence of transport, 
the dynamics of CBS can be used to extract microscopic 
parameters of the wave propagation in the disorder, in 
particular the scattering time and phase function [9HII]. 

The connection between CBS/WL, and strong, i.e., 
Anderson localization (AL) [12], which is exponential in 
real space, has, to our knowledge, remained experimen¬ 
tally unexplored, even though AL has also been observed 
in many systems nsuni. Theoretically, it was recently 
discovered that in Fourier space, besides the CBS peak 
at —ko, a twin peak at ko, called coherent forward scat¬ 
tering (CFS) peak [T5I1T8] should emerge at long times 
due to time reversal symmetry. For systems not spatially 
bounded, the symmetrization of the momentum distribu¬ 
tion with the CFS peak appearance has been proposed 
as a signature of strong AL [MIIH]. 

In this Letter, we explore the crossover between weak 
and strong localization, using a 2D photonic disorder gen¬ 
erated in a photorefractive crystal, a system where AL 
has already been studied nann]. We report the first 
observation of transverse CBS in a photonic disorder, at 
finite propagation time set by the crystal length, with a 
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FIG. 1: (color online) a : Sketch of the experiment. A quasi¬ 
plane wave with transverse momentum koy is launched in a 
^-invariant disorder and recorded at crystal output in Fourier 
space on a CCD camera, b : Principle of CBS : all pairs of 
scattering paths ^ —k^ (solid line) interfere constructively 
with their time-reversed counterpart (dashed line). 


peak contrast approaching the ideal value of 1. Further, 
varying in controlled manner the disorder strength Vfi 
and the incoming plane wave momentum ko, we study 
the full momentum distributions and the dynamics of 
the CBS peak, finding a non-monotonous evolution. In 
a first stage, the CBS peak grows, while the incoming 
wave peak decreases, leading to an asymmetric momen¬ 
tum distribution. This asymmetry reaches a maximum 
and for stronger disorder, it decays, as expected in the 
AL regime when the CFS peak appears cans]. Our data 
do not permit a definite observation of a CFS signal, but 
we show that the resymmetrization of the momentum 
distribution is compatible with the onset of AL in our 
system. Comparison to numerical calculations, showing 
similar features, supports this scenario. 

Our experiments. In our setup diiia, we realize a 
computer controlled, non-diffracting disorder using spa¬ 
tial light modulators (SLM), as explained in [20]. The 
spectrum of the writing beam is contained in a thin ring 







2 



B 


1 

b1 

I 

1 

d1 

I 

1 

1 

1 

d2 

I 

a3 t 

-Vr 

fi 

c3 1 

L.. 

d3 I 

a4 i 

1 1 

1 1 

1 1 

1 1 

1 1 

c4. 

1 1 

1 1 

d4 * 

1 1 

1 1 

1 1 

_1_1_ 

' ■ "y 1 ^ V “ w w 1 

1 1 

1 1 

_1_1_ 

LVaAw 

1 1 

1 1 

_1_1_ 

05. 

1 1 

_1_1 


FIG. 2: (color online) Observation of CBS peak formation 
(data A) for transverse momenta as indicated on columns, 
and disorder strengths (from first to last lines) VnlEa- = 
0,0.83,1.53,2.12,2.62. A : full pictures in far-held. Dashed 
circles show the elastic scattering circles. B : Azimuthal pro- 
hles. Dashed vertical lines show the positions of the incoming 
wave at 0i — j2 and the CBS peak at O 2 = 37r/2. 


(hltered with an SLM in Fourier space), ensuring that 
the disorder speckle held is almost ^-invariant through 
our 10-mm long SBN:75 crystal (see Fig. [^a), so that 
the photo-generated disordered refractive index Sn{x^ y) 
is also ^-invariant. The writing beam, at A = 27t jki, = 
532nm, is applied during variable times tw, which al¬ 


lows to adjust the disorder strength Vr uniiii], dehned 
as the variance of 6 n{x,y). To calibrate Vr, we use our 
method m m- During writing, we apply a bias held is 
Eq = +2kV/cm, so that the disorder is of focusing type 
(analogous to a red-detuned speckle potential for atoms). 
Due to the photorefactive effect [211423] . our disorder is 
anisotropic, with stronger scattering and localization ef¬ 
fects in the c-axis direction y. 

After writing, we send through the crystal a low in¬ 
tensity, SLM controlled, quasi-plane wave probe beam of 
adjustable transverse momentum k =/coy ESI. Its slowly 
varying amplitude T(x,^), in the paraxial regime, obeys 
a (2+l)D Schrodinger equation [24] 

= ( 1 ) 

oz 2po ^0 

where /3o = 27rno/A, is the transverse 

laplacian, z t plays the role of time t, and the potential 
V(x, y) is the refractive index : V(x, y) ^ —Sn{x, y) |34| . 

All data are averaged over ^ 40 disorder realizations. 
In this work we use a disorder with correlation length 
ctr = 6.53/im (same dehnition as for standard speckle dis¬ 
orders with a disk pupil function [23|26]), which dehnes 
a characteristic "time” la- = rieki^a^ = 1.17mm (along 
z), and an "energy" Ecr = {k^ne(T^)~^ = 7.2 x 10“^ [25] . 
which we use to normalize the disorder strength Vr. 

Observation of transverse CBS. Figurej^A shows far- 
held images of the probe beam after propagation in the 
transverse disorder, for various Vr and momenta ko. For 
increasing Vr, the initial peak at ko decreases, and one 
observes the formation, hrst, of a ring of elastic scatter¬ 
ing, then, of a very clear CBS peak at —ko, whose con¬ 
trast relative to the incoherent background in the scatter¬ 
ing circle approaches the ideal value of 1 in many pictures 
(see, e.g., a4, b4, a5-d5, and color bar). 

For quantitative analysis, we extract azimuthal prohles 
n{0), integrating the density on a ring of radius dk = 
O.bko (see Fig. [^a). Prohles n{0) are then htted with a 
two-peak function ES] 

n(^)=no(l+c,e-"W^+c,e-^), (2) 

where di = 7r/2 and 62 = 37r/2. In Fig. |^b (data c3 of 
Data B), the contrast of the CBS peak is C 2 = 0.9. It is 
important however to note that in our system, spectral 
broadening is strong, so that the scattering circle soon 
resembles a disk (see Fig. |^and[^. Thus, the choice of 
dk, which is arbitrary, inhuences the numerical outcome 
for n{0), and the ht. For example. Fig. [^c illustrates the 
dependence of C2 on dk (and on I/r). Nevertheless, we 
have checked that the choice of dk does not qualitatively 
change our further observations and conclusions. It can 
also be noted that the background no of incoherent scat¬ 
tering is not isotropic, as in Eq. In particular, for the 
larger ko one notes clear lobes of preferential scattering 
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directions (see Fig. |2}B .c3-d5), related to the particular 
autocorrelation function of our disorder [II 1 E 7 H 29 ]. 



FIG. 3: (color online) Construction of azimuthal profile n{0), 
for picture c3 in Data B. a : Full picture. The dashed area is 
the integration ring of radius dk = 0.5/co- b: Resulting n{0) 
with two-peak fit (smooth line), c : CBS peak contrast C 2 as 
function of dk/ko, for Data B, c2-c6 (/coctr = 0.98). 

Dynamics of momentum distribution asymmetry. To 
explore the crossover form weak to strong localization, 
we now focus on the asymmetry of the momentum dis¬ 
tribution. Several works have reported momentum dis¬ 
tributions where the asymmetry is opposite to the initial 
situation, due to the presence of a CBS peak at —ko on 
top of an isotropic scattering ring (e.g., 0130], and for 
simulations nni El). However, for strong localization 
(AL), one expects a resymmetrization of the momentum 
distribution, which should converge towards a twin peak 
structure, with CBS and CFS peaks of identical contrast 
and width m- 

In Fig. we show a second data set (B), where the 
probe beam has a 16 times wider area [33], and Vr 
reaches higher values. In a first stage, we observe the 
formation of well asymmetric profiles with a CBS peak 
dominating over the initial ko peak, and C 2 > ci (see, 
e.g., a4, a5, b3, b4, c4, c5). For larger Hr, the CBS peak 
decreases and the distribution recovers symmetry, with 
Cl C 2 . (see, e.g., a6, b6, c6). 

To quantify these observations. Figure gathers the 
evolutions of two indicators of asymmetry : the ratio of 
contrasts C 2 /C 1 , and the ratio n_/n+ of the total den¬ 
sity in the lower {ky < 0) vs upper half {ky > 0) of the 
momentum space. The first ratio indicates the relative 
weight of —ko vs ko components, independently from a 
decrease of both ci and C 2 at strong disorder (as seen in 
Fig. ic and Fig. [^. The second ratio is more general. 
To compare all our data at different ko, we rescale the 
Vr using the value that we define as the value for 
which the amplitude ni = noCi of the initial peak has 
decayed to 10% of its initial value. This rescaling makes 
all curves for different ko collapse onto a typical path. 
Physically, this procedure removes the /c-dependence of 
the scattering mean free paths ls{k), i.e., one essentially 
has V^{k) /s(^)~^- In data A (Fig. |^a, d), both asym- 


A ko‘^R= 0-43 k„°R= 0-68 kQa^= 0.98 k„a^= 1.38 
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FIG. 4: (color online) Dynamics of CBS formation (data 
B) for disorder strengths (from first to last rows) Hr/Fct = 
0,1.53,2.12,2.62,3.04,3.4. 
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Data A 




Data B Simulations 






Simulations 



FIG. 5: (color online) Dynamics of the asymmetry of contrasts C 2 /C 1 (a-c) and of the momentum distribution n+/n_ 
(d-f) in data sets A, B, and simulations, for disorder strengths Vr normalized to Vr (see text). Data A : Vr/Eo- = 
1.41,1.63,1.76,1.84,1.94,2.12,2.45. Data B : Vr/E^ = 1.4,1.35,1.7,2.15. Simulations : Vr/E^ = 0.51,0.61,0.74,0.85. (g-1) : 
Examples of simulations results for koaR = 1. 


metry ratios increase. In data B (Fig. [^b, e), both reach 
a maximal value (about 2 for C 2 /C 1 and 1.07 for n+/n_), 
and then decrease towards approximately 1, meaning full 
resymmetrization at the maximal Vr. 

Comparison to numerical calculations. To test our 
observations, we simulate Eq. using a split-step beam 
propagation algorithm, and a gaussian probe beam of 
r.m.s width Axsim = 200/im. For the disordered 6n{x, y)^ 
we solve the standard steady-state anisotropic equations 
for the photorefractive effect (see m or Eq. 9 in |22)). 
Yet, an important difference should be noted. In the ex¬ 
periments, the disorder strength Vr is varied by adjusting 
the photorefractive writing time, while the writing inten¬ 
sity Jw is constant and much larger than the saturation 
intensity (/wZ-^sat r\j 10^) [21]. Whereas, in the numerics, 
steady-state is assumed (tw = 00) and Vr is adjusted by 
varying Jw/^sat = 0.1 — 5, i.e., from weak to considerable 
conditions of photorefractive saturation. 

We average the simulations over 50 disorder realiza¬ 
tions, and perform the same data analysis as for the ex¬ 
periments. Overall, the simulations (Fig. |^c,f) display 
similar trends as the experiments for both asymmetry 
ratios, with first a maximum (also ^ 2 for C 2 /C 1 , but 
a lower ^ 1.02 for n+/n_), and then a decay, which is 
much slower than in the experiment. In Fig. [^g-1, we 
show examples of the simulation results for koaR = 1. 
In the first case (g,j), the ko peak dominates. In the 
intermediate case (h,k), the CBS peak dominates, and 
maximal asymmetry is reached. In the third case (i,l), 
the asymmetry has almost disappeared. There, one also 
notes a strong disorder broadening, as also seen in the 
experiment (see Fig. iA. c6), and a strong anisotropy 


in the momentum distribution, more pronouced than in 
the experiment. Despite several quantitative differences, 
which may well be due to the above-mentioned difference 
of disorder type and saturation effects, one notes that the 
maximum asymmetry occurs for similar absolute values 
of Vr, for example, for C 2 /C 1 and ko — 1, at Vr 2.4E'cr 
in data B and Vr ^ 2.6Efj in simulations. In both cases 
also, the maximum of n^ju- occurs for slightly larger 
Vr. And we neither found any sharp CFS peak in the 
simulations. Most importantly, the similarity between 
measured and simulated dynamics, allows us to consider 
the resymmetrization as a fully coherent effect, given that 
lattice imperfections (for example residual z-variations of 
the disorder), are absent in the numerics. 

Can we interpret the resymmetrization at stronger dis¬ 
order to the onset of strong (Anderson) localization ? 
One could firstly note that exponential localization in real 
space, interpreted as AL, has been well observed in essen¬ 
tially the same system for similar parameters [1^1 HOI [Hj- 
Analysing our parameters more quantitatively brings fur¬ 
ther insight. Our measurements are all carried at fixed 
propagation ’Time" t /ter — 8.55. In the numerical work 
m, using a similar disorder (although the pupil func¬ 
tion is a disk instead of a ring in our case), the timescale 
for the CFS peak contrast to reach 0.2 is t/ta- ^ 20 for 
Vr = 3 — 4F^cr, and koaR = 1.5. Following our rescaling 
procedure, in our Data B with similar Vr = 3.4crR, for 
/coctr 0.8, one expects similar effects for t/to- ^ 11 [36] . 
which is fairly compatible with our parameters. It thus 
seems very reasonable to attribute the resymmetrization 
of the momentum distribution at finite time in our data, 
to the crossover form weak to strong localization. 
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Conclusion. We have reported the first observation 
of transverse CBS and full 2D momentum distributions 
in a photonic disorder. We studied the dynamics of the 
CBS peak and of the momentum distribution asymme¬ 
try. We noted, first, a growth of the CBS peak, whose 
contrast approaches 1, and then for stronger disorder, a 
resymmetrization, which we attributed to the onset of 
Anderson localization. Further analysis could provide 
better understanding of scattering, spectral broadening, 
and coherence in our disorder, and better estimates of 
the localization parameters, which may allow to resolve 
the CFS peak appearance in our system. 
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